Frequency comparisons and absolute frequency measurements of 
171 Yb + single-ion optical frequency standards 

E. Peik 1 '*, B. Lipphardt 1 , H. Schnatz 1 , T. Schneider 1 , Chr. Tamm 1 , 

S. G. Karshenboim 2,3 

1 Physikalisch-Technische Bundesanstalt, 
Bundesallee 100, 38116 Braunschweig, Germany 
2 D. I. Mendeleev Institute for Metrology (VNIIM), 198005 St. Petersburg, 

Russia 

3 Max-Planck-Institut fiir Quantenoptik, 85748 Garching, Germany 
* e-mail: ekkehard.peik@ptb.de 

Abstract 

We describe experiments with an optical frequency standard based on 
a laser cooled m Yb + ion confined in a radiofrequency Paul trap. The 
electric-quadrupole transition from the 2 Si/ 2 (F = 0) ground state to the 
2 Ds,/2{F = 2) state at the wavelength of 436 nm is used as the reference tran- 
sition. In order to compare two 171 Yb + standards, separate frequency servo 
systems are employed to stabilize two probe laser frequencies to the reference 
transition line centers of two independently stored ions. The experimental 
results indicate a relative instability (Allan standard deviation) of the optical 
frequency difference between the two systems of <7 y (1000s) = 5 • 10~ 16 only, 
so that shifts in the sub-hertz range can be resolved. Shifts of several hertz 
are observed if a stationary electric field gradient is superimposed on the ra- 
diofrequency trap field. The absolute optical transition frequency of Yb + at 
688 THz was measured with a cesium atomic clock at two times separated by 
2.8 years. A temporal variation of this frequency can be excluded within a 
la relative uncertainty of 4.4 • 10 -15 yr" 1 . Combined with recently published 
values for the constancy of other transition frequencies this measurement 
provides a limit on the present variability of the fine structure constant a at 
the level of 2.0 ■ 10~ 15 yr" 1 . 
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1 Introduction 



Time and frequency are presently by far the most precisely measurable physi- 
cal quantities and optical atomic frequency standards hold great potential for 
further improvements in stability and accuracy. A single trapped and laser 
cooled ion can be used as a reference that offers excellent control over possible 
systematic effects. Presently, a number of groups is pursuing research along 
these lines with different ions (see PQ E] for recent reviews). High- resolution 
spectroscopy and precise frequency measurements using femtosecond laser 
frequency comb generators [3] have been performed with reference transi- 
tions in Sr+ H|3], In + 013, Yb+ CHI HI] and Hg+ |I2l HSJ 

In this paper we review two recent experiments at PTB with an optical 
frequency standard based on the m Yb + ion: (i) frequency comparisons be- 
tween two ions in separate traps are used to investigate the stability and to 
evaluate systematic frequency shifts [Hj, (ii) a sequence of absolute optical 
frequency measurements has been made using a cesium fountain clock in or- 
der to test the reproducibility of the standard and to obtain a stringent limit 
on a possible temporal variation of the fine structure constant [To] . 

2 171 Yb + optical frequency standard 

The 171 Yb + ion is attractive for an optical frequency standard because ref- 
erence transitions with vanishing low-field linear Zeeman frequency shift 
are available in a level system with relatively simple hyperfine and mag- 
netic sublevel structure (cf. Fig. la). The electric quadrupole transition 
( 2 Si/2,F = 0) — > ( 2 D 3 /2,F = 2) at 436 nm wavelength can serve as the 
reference transition, with a natural linewidth of 3.1 Hz. A single ytter- 
bium ion is trapped in a miniature Paul trap and is laser-cooled to a sub- 
millikelvin temperature by exciting the low-frequency wing of the quasi-cyclic 
(F = 1) — > (F = 0) component of the 2 Si/ 2 — > 2 Pi/2 resonance transition at 
370 nm. A static magnetic field of approximately 0.3 mT is applied in order 
to prevent optical pumping to a nonabsorbing superposition of the magnetic 
sublevels of the F = 1 ground state. The natural linewidth of this transition 
is 21 MHz, which implies a one-dimensional kinetic temperature of 0.6 mK at 
the Doppler cooling limit. A weak sideband of the cooling radiation provides 
hyperfine repumping from the F = ground state to the 2 P\/2{F = 1) level. 
At the end of each cooling phase, the hyperfine repumping is switched off in 
order to prepare the ion in the F = ground state. The rapid spontaneous 
decay from the 2 P\/2 state to the (F = 2) sublevel of the metastable 2 D 3 / 2 
state that occurs during laser cooling is compensated for by coupling this 
level with an additional laser at 935 nm wavelength to the [3/2] 1/2 (F = 1) 
state, from where the ion readily returns to the ground state. The (F = 2) 
sublevel of the 2 D 3 / 2 state is not rapidly populated or depleted by the laser 



2 



cooling excitation. Individual quantum jumps to this state due to excita- 
tion of the reference transition can therefore be detected using the electron 
shelving scheme. 




Figure 1: a) Simplified level scheme of the 171 Yb + ion. b) Excitation spec- 
trum of the 2 Si/ 2 (F = 0,m F — 0) — > 2 D 3 / 2 (F = 2,mp = 0) transition, 
obtained with 20 probe cycles for each value of the detuning. The linewidth 
of 10 Hz is approximately at the Fourier limit for the employed probe pulses 
of 90 ms duration. 

The S — D electric quadrupole reference transition is probed by the fre- 
quency doubled radiation from a diode laser emitting at 871 nm. The short 
term frequency stability of this laser is derived from a temperature-stabilized 
and seismically isolated high-finesse reference cavity. Figure lb shows a high- 
resolution excitation spectrum obtained with 90 ms long laser pulses, leading 
to an approximately Fourier-limited linewidth in the range of 10 Hz, or a reso- 
lution Au/u of 1.4- 10~ 14 . Since the duration of the probe pulse is longer than 
the lifetime of the excited state (51 ms), the observed maximum excitation 
probability is limited by spontaneous decay. In order to operate the system 
as a frequency standard, both wings of the resonance are probed alternately, 
and the probe light frequency is stabilized to the line center according to the 
difference of the measured excitation probabilities. 



3 Frequency comparisons between two trapped 
ions 

To enable a quantitative study of systematic frequency shifts on the sub-hertz 
scale, we compare the line center frequencies of two 171 Yb + ions stored in sep- 
arate traps. The scheme of the frequency comparison experiment is shown in 
Fig. 2. Both traps use the same cooling laser setup and synchronous timing 
schemes for cooling, state preparation, and state detection. The beam from 
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Figure 2: Experimental setup for the comparison of two 171 Yb + frequency 
standards. AOM: acusto-optic modulators, providing independent frequency 
shifts between the probe laser and the two ion traps. 

the probe laser is split and two separate frequency shift and servo systems are 
employed to stabilize the probe frequencies to the reference transition line 
centers of the two ions. The atomic resonance signal is probed in typically 
eight excitation cycles until an error signal is calculated and the detunings 
between the probe laser frequency and the probe light beams incident on the 
traps are corrected. The servo time constants are in the range of 10 s. In 
order to minimize servo errors due to the drift of the probe laser frequency, 
a second-order integrating servo algorithm is used. The differences of the 
detunings imparted on the probe beams are averaged over time intervals of 
1 s and recorded. As a measure of the instability, the Allan deviation of a 
typical data set is shown in Fig. 3. The result indicates a relative instability 
of the frequency difference between the two systems of o" y (1000 s) = 5 ■ 1CT 16 . 
The atomic resonance signals were resolved with Fourier-limited linewidths 
of approximately 30 Hz in both traps. For averaging times that are longer 
than the time constants of the servo systems, the so-called quantum projec- 
tion noise JH] limits the instability. This contribution to the instability is 
dominant because the state measurement is done on a single atom only. We 
have performed a numerical calculation which simulates the effect of quan- 
tum projection noise for the realized experimental conditions and include the 
result as the solid line in Fig. 3, showing good agreement with the measured 
data. 

The mean frequency difference between the two traps was generally in the 
hertz or sub-hertz range for nominally unperturbed conditions [T3|. The 
dominant source of systematic uncertainty in the 171 Yb + optical frequency 
standard is the so-called quadrupole shift of the atomic transition frequency. 
It is due to the interaction of the electric quadrupole moment of the D 3 / 2 
state with the gradient of the static electric field in the trap. In order to 
experimentally determine the 171 Yb + quadrupole moment, a static field gra- 
dient was generated in one of the traps by superimposing a constant voltage 
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Figure 3: Allan standard deviation of the frequency differences observed 
between the two traps, normalized to the optical frequency of 688 THz. The 
dashed line connects the experimental points and the solid line shows the 
result of a Monte Carlo simulation of the servo action for the case that 
the fluctuations of the atomic resonance signals are determined by quantum 
projection noise. 

on the radiofrequency trap drive voltage. The orientation of this field gra- 
dient is determined by the symmetry axis of the trap. The other trap was 
operated with a pure rf voltage and served as a reference. From this exper- 
iment, we infer a quadrupole moment of Q = (3.9 ± 1.9)eag for the 2 D^/2 
level with e being the elementary charge and ao the Bohr radius. The un- 
certainty is mainly determined by the uncertainty of the measurement of the 
angle between the magnetic field and the trap axis. Depending on the orien- 
tation, a field gradient of 1 V/mm 2 may produce a frequency shift of up to 
6 Hz. Averaging the frequency over three orthogonal quantization axes - i.e. 
axes of the magnetic field - allows to determine the unperturbed transition 
frequency [T7| . 

4 Absolute frequency measurements 

The absolute frequency of the m Yb + standard at 688 THz was measured 
relative to a cesium fountain clock ^H] > using the setup shown schematically 
in Fig. 4a. A hydrogen maser is referenced to the cesium clock and delivers a 
signal at 100 MHz with a relative instability below 2 ■ 10~ 15 in 1000 s of aver- 
aging. The link between the optical and microwave frequencies is established 
with a femtosecond-laser frequency comb generator. The m Yb + frequency 
was measured first in December 2000 and January 2001. The result was 
688 358 979 309 312 Hz with a total la measurement uncertainty of 6 Hz [8]. 
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Figure 4: a) Schematic of the setup for the absolute frequency measurements 
of the 171 Yb + frequency standard, b) Results of absolute frequency measure- 
ments of the m Yb + standard versus the date of the measurement (MJD: 
Modified Julian Date). 

In November 2003, the frequency was remeasured after improvements in the 
resolution of the ionic resonance and in the frequency comb generator setup. 
The new value has a total la uncertainty of 6.2 Hz and is insignificantly lower 
by 2.3 Hz. The dominant source of systematic uncertainty is the quadrupole 
shift that may be produced by patch charges on the trap electrodes. The 
static quadrupole potential associated with these charges will also alter the 
oscillation frequencies of the ion in the trap. Measurements of the oscillation 
frequencies allow to estimate the uncertainty due to the quadrupole shift as 
3 Hz under the conditions of the experiment. The total uncertainty of 6.2 Hz 
(corresponding to a relative uncertainty of 9 ■ 10~ 15 ) contains also a system- 
atic contribution of 2.2 Hz from the microwave reference and frequency comb 
generator and a total statistical contribution of 4.8 Hz. 

From the sequence of absolute frequency measurements (cf. Fig. 4b) we 
deduce a value for the fractional temporal variation of the frequency ratio 
fyb/fcs °f ( — 1-2 ±4.4) ■ 10~ 15 yr _1 , consistent with zero. This result can 
be used to obtain a limit on the temporal variation of the fine structure 
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constant. 



5 Variations of fundamental constants 

Discussions about a possible temporal variability of fundamental constants 
of nature date back until the 1930s, with Dirac's large number hypothesis 
[T9"j being a prominent example. Meanwhile, strong motivation to search for 
such an effect has developed because theories that attempt to unify the fun- 
damental interactions allow or even imply temporal variations of the coupling 
constants (see e.g. [20 ). Variations of the constants have been searched for in 
various systems [2UI22]- Sommerfeld's fine-structure constant a = e 2 /47re frc 
(composed of the elementary charge e, Planck's constant h, the speed of light 
c, and the electric constant of vacuum eo) is presently the most important 
test case because it is the fundamental coupling constant of the electro- 
magnetic interaction. As a dimensionless quantity a ~ 1/137.036 can be 
determined without reference to a specific system of units. This point is im- 
portant because the realization of the units themselves may also be affected 
by variations of the constants. 

An indication for a possible variation of a comes from a shift of wavelengths 
of specific absorption lines produced by interstellar clouds in the light from 
distant quasars These observations seem to suggest that about 1 ■ 10 10 yr 
ago the value of a was smaller than today by Aa/a = (— 0.543±0.116) • 10~ 5 , 
representing 4.7a evidence for a varying a |24j . Assuming a linear increase in 
a with time, this would correspond to a drift rate dlna/dt = (6.40 ± 1.35) ■ 
10~ 16 yr -1 [23] • Recent evaluations of new data on quasar absorption lines 
are consistent with Act = for all look-back times [2*5] . 

To perform a laboratory tests for the variability of a in the present epoch 
one can compare different electronic transition frequencies in the optical 
range j2Hl I2*7j . For the analysis, we express the electronic transition fre- 
quency as 



where Ry = m e e 4 /(8eo/i 3 ) — 3.2898 ■ 10 15 Hz is the Rydberg constant in 
SI units of frequency, appearing here as the common atomic scaling factor. 
C is a numerical constant which depends only on the quantum numbers 
characterising the atomic state and which is independent of time. F(a) 
is a dimensionless function of a that takes into account level shifts due to 
relativistic effects |2H] . In many-electron atoms these effects lead to different 
sensitivities of the transition frequencies to a change of a, with the general 
trend to find a stronger dependence in heavier atoms [29] . The relative 
temporal derivative of the frequency / can be written as: 




(1) 
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dt ~ dt 
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The first term din Ry/dt represents a variation that would be common to all 
measurements of electronic transition frequencies: a change of the numerical 
value of the Rydberg constant in SI units, i.e. with respect to the 133 Cs hy- 
perfine splitting. We distinguish here between the numerical value of Ry and 
the physical parameter m e e 4 //i 2 which determines the atomic energy scale. 
A change of the latter would not be detectable in this type of measurement 
because the 133 Cs hyperfine splitting also scales with this quantity. Conse- 
quently, the detection of a non-zero value of din Ry/dt would have to be 
interpreted as being related to a change of the nuclear magnetic moment of 
133 Cs. The second term in Eq. 2 is specific to the atomic transition under 
study. The sensitivity factor A for small changes of a has been calculated 
for a number of cases of interest by Dzuba, Flambaum et al. [2H]- Assuming 
that all possible drifts are linear in time over the duration of the frequency 
comparisons, the fractional drift rate dlna/dt can be obtained from at least 
two measured drift rates of transition frequencies din f/dt if these have been 
measured against a common reference, and if the sensitivity factors A for the 
two transitions are different. 

To derive the limit on dlna/dt we combine the data from 171 Yb + with 
a published limit on the drift rate of a transition frequency in 199 Hg + [T3*] . 
obtained at NIST (Boulder, USA). In Hg + the transition 5<i 10 6s 2 Si/2 — > 
5d 9 6s 2 2 D 5 /2 at 282 nm (1065 THz) with a natural linewidth of 1.9 Hz serves 
as the reference. A first absolute frequency measurement of this transition 
with respect to a cesium clock was published in 2001 with an uncertainty of 
9 • 10 -15 [30 . A sequence of measurements over the period from August 2000 
to November 2002 has resulted in a constraint for the fractional variation of 
/kg/ics at the level of 7 ■ 10~ 15 yr _1 (lcr) 13 . The sensitivities of the ytter- 
bium and mercury transition frequencies to changes of a are quite different: 
A Yh = 0.88 and A Ug = -3.19 [2E|. 

A further precision test for the constancy of an optical transition frequency 
was presented by Fischer et al. who measured the IS* — > 2S transition 
frequency in atomic hydrogen in two periods in July 1999 and in February 
2003, using a transportable cesium fountain frequency standard jSIj. The 
limit on the fractional drift rate for the hydrogen frequency that is deduced 
from these experiments is (—3.2 ± 6.3) • 10 -15 yr _1 with the sensitivity 
factor An being zero. 

Plotting the three limits on the frequency drift rates as a function of A 
(Fig. 5a) and performing a weighted linear regression we obtain from the 
slope of the fit: 

din a , i 

^ r = (-0.3±2.0)-HT 15 yr- 1 (3) 

for the present value of the temporal derivative of the fine structure constant 
at the confidence level of lcr. This limit is based on a model-independent 
analysis, not using assumptions on correlations between drift rates of different 
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Figure 5: a) Relative frequency drift rate versus sensitivity factor A for the 
S — > £) reference transitions in the Hg + and Yb + ions and for the IS — > 2S 
transition in atomic hydrogen. The straight line is the result of a weighted 
linear regression. A significant deviation of the slope from zero would indicate 
a non-zero time derivative of the fine structure constant a. b) Constraints 
on the variability of a and the numerical value of Ry as obtained from the 
experimental points in (a). Stripes mark the la uncertainty of the individ- 
ual measurements and the central hatched area is bounded by the standard 
uncertainty ellipse resulting from the combination of all three experiments. 

quantities and not postulating the drift rate of any quantity to be zero. To 
our knowledge, this is so far the most stringent limit on da/dt obtained from 
atomic clock comparisons [15] . Another related analysis based on the data 
from hydrogen and mercury is presented in [31j . 

The intercept of the straight line representing the dependence of the optical 
frequency drift rate on A with the line A = in Fig. 5a can be used to 
determine the drift rate of the numerical value of the Rydberg frequency. 
This quantity is of great metrological importance. A non-zero value would 
imply among other consequences that a time scale generated with a cesium 
clock would not be uniform with respect to other atomic oscillations. The 
result that we obtain here is consistent with zero: 
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Fig. 5b shows how the measurements of the three transition frequencies 
contribute to the constraints on the temporal derivatives of a and Ry. The 
central hatched area is bounded by the standard uncertainty ellipse (see Ref. 
[THj for details). 

With these results we reach a sensitivity for Aa/At with At ~ 3 yr that 
approaches that of the quasar observations [2U 122] to within about a factor 
of three if a linear change of a is assumed on a cosmological time scale. 
Further improvements in the data from clock comparisons can be expected 
since trapped ions have the potential to reach relative uncertainties in the 
range of 10~ 18 for the comparison of optical transition frequencies pQ. The 
kind of analysis that we have applied here to the frequencies of the S —> D 
quadrupole transitions in 171 Yb + and 199 Hg + and of the IS — > 2S transition 
in H may be extended to other atomic systems as soon as precision data on 
frequency drift rates become available. 
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